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Abstract 

For physical systems described by smooth, finite-range and confining microscopic interac- 
tion potentials V with continuously varying coordinates, we announce and outline the proof 
of a theorem that establishes that unless the equipotential hypersurfaces of configuration space 
Ti v = {(<7i, • • • , <?tv) G M. N \ V(qi, . . . , qN) = v}, v £ R, change topology at some v c in a given inter- 
val [i>o,i>i] of values v of V, the Helmoltz free energy must be at least twice differentiable in the 
corresponding interval of inverse temperature (/3{vq), P(v\)) also in the N — > oo limit. Thus the 
occurrence of a phase transition at some (3 C = (3(v c ) is necessarily the consequence of the loss of 
diffeomorphicity among the {T, v } v<Vc and the {T lv } v>Vc , which is the consequence of the existence 
of critical points of V on S-y=v c j that is points where VV^ — 0. 
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Phase transitions (PTs) are phenomena which bring about qualitative physical changes 
at the macroscopic level in presence of the same microscopic forces acting among the con- 
stituents of a system. Their mathematical description requires to translate into quantitative 
terms the mentioned qualitative changes. The standard way of doing this is to consider 
how the values of thermodynamic observables, obtained in laboratory experiments, vary 
with temperature, or volume, or an external field, and then to associate the experimentally 
observed discontinuities at a PT to the appearance of some kind of singularity entailing a 
loss of analyticity. Despite the smoothness of the statistical measures, after the Yang-Lee 
theorem [l| we know that in the iV — > oo limit non-analytic behaviors of thermodynamic 
functions are possible whenever the analyticity radius in the complex fugacity plane shrinks 
to zero, because this entails the loss of uniform convergence in N (number of degrees of 
freedom) of any sequence of real-valued thermodynamic functions, and all this depends on 
the distribution of the zeros of the grand canonical partition function. Also the other devel- 
opments of the rigorous theory of PTs, like that due to Dobrushin, Lanford and Ruelle on 
Gibbs measures j2ly|, identify PTs with the loss of analyticity. 

However, we can wonder whether this is the ultimate level of mathematical understanding 
of PT phenomena, or if some reduction to a more basic level is possible. The present 
paper addresses just this point and aims at providing a non-technical presentation of a 
new rigorous result, reported in Ref. j^], making its conceptual meaning and prospective 
physical interest accessible without going through the details of a lengthy mathematical 
proof. The new theorem says that non-analyticity is the "shadow" of a more fundamental 
phenomenon occurring in configuration space: a topology change \^ within the family of 
equipotential hypersurfaces = {(qi,...,q N ) e M | V(q%, . . . , g^r) — v }, where V and 
qi are the microscopic interaction potential and coordinates respectively. This topological 
approach to PTs stems from the numerical study of the Hamiltonian dynamical counterpart 
of phase transitions, and precisely from the observation of discontinuous or cuspy patterns 
displayed by the largest Lyapunov exponent at the transition energy |f| (or temperature). 
Lyapunov exponents measure the strength of dynamical chaos and cannot be measured in 
laboratory experiments, at variance with thermodynamic observables, thus, being genuine 
dynamical observables they are only measurable in numerical simulations of the microscopic 
dynamics. To get a hold of the reason why the largest Lyapunov exponent Ai should probe 
configuration space topology, let us first remember that for standard Hamiltonian systems, 
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described by H = YL%=i §.Pi +^(<Zij ■ ■ ■ 5 Qn)j M is computed by solving the tangent dynamics 
equation 

^Hwj,/ =0 ' (1) 

where q(t) = ..,<?*(*)], and then A x = lim^ l/2t log(Sf =1 [C 2 (t) + ^(W^it^O) + 

^ 2 (0)]). If there are critical points of V in configuration space, that is points q c = . . . , g^] 
such that W(g)| = 0, according to the Morse Lemma |8j], in the neighborhood of any 
critical point q c there always exists a coordinate system q(t) = [qi(t), ..,gjv(t)] for which 

^(g) = V(q c ) -ql q\ + q 2 k+1 + ••• + &, (2) 

where k is the index of the critical point, i.e. the number of negative eigenvalues of the 
Hessian of V. In the neighborhood of a critical point, Eq.(J2J) yields df^V = ±5^ which, 
substituted into Eq.Q, gives k unstable directions which contribute 9] to the exponential 
growth of the norm of the tangent vector £. This means that the strength of dynamical 
chaos, measured by the largest Lyapunov exponent Ai, is affected by the existence of critical 
points of V. In particular, let us consider the possibility of a sudden variation, with the 
potential energy v, of the number of critical points (or of their indexes) in configuration 
space at some value v c , it is then reasonable to expect that the pattern of \\(v) - as well 
as that of Xi(E) since v = v(E) - will be consequently affected, thus displaying jumps or 
cusps or other "singular" patterns at v c (this heuristic argument has been given evidence in 
the case of the XY-mean-field model, see jsj and 0]). On the other hand, Morse theory |s| 
teaches us that the existence of critical points of V is associated with topology changes of the 
hypersurfaces {£„}„ e R, provided that V is a good Morse function (that is: bounded below, 
with no vanishing eigenvalues of its Hessian matrix). Thus the existence of critical points 
of the potential V makes possible a conceptual link between dynamics and configuration 
space topology, which, on the basis of both direct and indirect evidence for a few particular 
models, has been formulated as a topological hypothesis about the relevance of topology 
for PTs phenomena. In what follows, we show that, for a large class of physically meaningful 
potentials, this conjectural status of the art turns into a qualitatively new one because we 
can prove the following 

Theorem. Let Vjv(gi, . . . , qw) '■ ^ N — > K, be a smooth, bounded from below, finite-range 
and confining potential L r lC]. Denote by E„ := V^ x {v), its level sets, or equipotential 
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hypersurfaces, in configuration space. Then let v = v/N be the potential energy per degree 
of freedom. 

If there exists N , and if for any pair of values v and v' belonging to a given interval 
h = [vo, Vi] and for any N > N 

T, m is diffeomorphic to T>Nv' 

then the sequence of the Helmoltz free energies {F N ({3)} N£ ^ - where (3 = 1/T (T is 
the temperature) and (5 E Ip — ((3(vq), (3{vx)) - is uniformly convergent at least in C 2 (Ip) 
[the space of twice differentiable functions in the interval Ip], so that lirriAr^oo-F/v € C 2 (Ip) 
and neither first nor second order phase transitions can occur in the (inverse) temperature 
interval (j3(v ), (3(vi)). 

Where the inverse temperature is defined as (3(v) = dS^ N \v) / dv and S^ N \v) = 
N^ 1 log J v ^ <JjN d N q is one of the possible definitions of the microcanonical configurational 
entropy. The intensive variable v has been introduced to ease the comparison between 
quantities computed at different iV-values. 

This theorem means that a topology change of the {£„},,eR at some v c is a necessary condition 
for a phase transition to take place at the corresponding energy or temperature value. The 
topology changes implied here are those described within the framework of Morse theory 



through attachment of handles 



Remark 1. The topological condition of diffeomorphicity among all the hypersurfaces T,nv 
mth ~ v £ M Has an ana^ca, c_e: the absence of c M points o, V in the 
interval [vo,vi . This is proved in Lemma 1 of Ref\A\ by adapting to the £„ Bott's "critical 
neck theorem" which applies to the manifolds M v = {(qi, ^at) G M. N \V(qi, q^) < v}. 
Apart from this initial link with topology, the proof proceeds in the domain of Analysis. 

Remark 2. In the proof we resort to the concept of uniform convergence - from elementary 
functional analysis - of a sequence of functions, and to the fact that the limit of a sequence 
of smooth functions can be non-smooth. This way of tackling the thermodynamic limit is in 
the spirit of the celebrated Yang-Lee theorem 

Let us now outline the proof by focusing on the main ideas (details can be found in 
Under the crucial hypothesis of diffeomorphicity of the hypersurfaces f° r u 

[■Uo,Vi], we want to prove that the thermodynamic limit of the Helmoltz free energy, 
Foo(fi) = hniTv^oo Fn{/3), is at least twice differentiable, so that first or second order 
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phase transitions are absent. For standard Hamiltonians, each function Fn(/3) reads as 
Fn(P) = — (2/3)" 1 \og{ir/0) — /jv(/3)//3, sum of a part coming from the kinetic energy- 
term, and a configurational part /jv(/3) = (V-W)l°g/ ^9 ex P[ — /^(g)]. Thus, in order 
to prove that F^fS) G C 2 (/g), we have to show that the sequence of smooth functions 
{-F/v(/3)}jv e N is uniformly convergent at least in C 2 (Ip) in the limit iV — > oo, or equivalently, 
since {2f3)~ l log(7r//3) remains always smooth in the limit iV — > oo, we have to show that the 
sequence of smooth functions {/;v(/3)}./V€N is uniformly convergent in C 2 (Ip) when N — > oo. 
Now, at any AT, /jv(/3) is related to the microcanonical entropy Sj^ through the Legendre 
transform: Sj^ (?;) = Jn{P) + (3v. Hence, to prove that foo(P) is twice differentiable with 
respect to (3 we need to prove that Soo\v) is three times differentiable with respect to v. 

Eventually, we consider the equivalent definition, at large N, of the configurational mi- 
crocanonical entropy 

S N (v) = hog n{v,N) = ±=]og [ J^,, (3) 



N N \ VV 

** L-iV II I i 

which also implicitly defines Q as the surface integral in the r.h.s., where da is the (N — 1)- 
dimensional surface element of E„, and where ||W|| = [^iLiidV/dqi) 2 ] 1 / 2 ; since Sn(v) 
has the same thermodynamic limit of the entropy S^ N (v), that is sL\v) = Soo(v), we are 
left with the problem of proving that the sequence of smooth functions {Sn(v)}n<en [where 
Sn(v) — Sn(vN)}, is uniformly convergent in C 3 ^), the space of three times differentiable 
functions in the interval J s , in the limit iV — > oo. The reason for using Sn{v) instead of 



Sfc\v) will be soon clear. After the Ascoli theorem 



12j, in order to prove that S^v) is 



three times differentiable, we need to prove that for v G 1% — \vq,Vi] and for any N, the 
function Sn{v) and its first four derivatives are uniformly bounded in N from above, that 
is, for any iV G N and v G [vq, Vi) 

d k S N 



sup \Sn(v)\ < oo , sup 



< oo , k = 1, ..,4. (4) 



dv k 

We prove the Theorem by proving that these bounds are the consequence of the diffeo- 
morphicity among the T*m, for v G [v , Vi). 

From Eq.(j2I) the first four derivatives of S^(v) are trivially computed: 

dS N _ 1 Q'{v,N) dv _ Q'(v,N) 

dv ^ ~ N Q(v,N) ' dv ~ n(v,N) ^ ' 

and, using a compact notation, d$S N = N[tt"/Q- (fi'/fi) 2 ], SgSW = N 2 [Vl"' /fi-3fi"fi' /ft 2 + 
2(Q'/n) 3 } &nddlS N = N 3 [n iv /n-4n m Q , /Q 2 -3(Q''/Q) 2 + 12n''(n , ) 2 /n 3 -6(Q'/Q)% where 



the prime indexes stand for derivations of Q(v, N) with respect to v — vN. In order to verify 
whether the conditions (jU are fulfilled, we must be able to estimate the iV-dependence of 
all the addenda in these expressions for the derivatives of Sjy. 



Being the assumption of diffeomorphicity of the £jvd equiva. 



points of the potential, we can use the derivation formula 
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ent to the absence of critical 



= / IIWII A' 



1 



w 



da 



W 



(6) 



where A k stands for k iterations of the operator 



A{.) 



W_ 
IWI 



l|W|| 



The technical reason to work with Sn instead of is now evident: the derivatives of 
fl(v, N) are transformed into the surface integrals of explicitly computable combinations 
and powers of a few basic ingredients, like ||W||, dV/dqi, d 2 V/dqidqj, d 3 V/dqidqjdqk and 
so on. This is a technically crucial step to prove the Theorem. 

The first uniform bound in Eq.(@J, |SV(^)| < oo, is a simple consequence of the intensivity 
of S N (v). 

To prove the boundedness of the first derivative of Sn, we first compute its expression 
by means of Eqs.(jH} and ©, which reads 



OS 



N 



dv 



1 



AV Jlij&VdlV&iV 



da 



NV 



with diV = dV/dq l and i, j 
dS N 



dv 



< 



W|| 2 ||vy|| 4 

, N, whence (with an obvious meaning of (-)t, v 



(7) 



+ 2 



1j v 



IWI 



(8) 



The A^-dependences of the derivatives of Sn are estimated at constant potential energy 
density v, for any given v G [v~o, Vi], thus we can think of increasing by glueing together 
an increasing number k of replicas of a given building block of N particles at potential energy 
v = vN . Now the hypothesis of diffeomorphicity of the plays again a crucial role, in fact 
in the absence of critical points, for each building block we have min || Wat || 2 > C 2 > 0, 
with C a constant, and as only short-range potentials are considered, the larger iVo and 
N(= kNo) the less relevant are the boundary interactions among the building blocks. Thus, 
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at large N, min || W^|| 2 > C 2 N, where C = C/Nq is a constant; for an upper bound estimate 
of Eq.fJSJ we replace in its denominators the lower bound C 2 N of min ||W^|| 2 



OS 



N 



dv 



< 



<|AV|> 



+ 2 



C 2 iV C 4 N 2 

where now we have to estimate the ^-dependence of the numerators. To this purpose, as we 
have assumed that V is smooth and bounded below, we note that (| AV \)x v = (\ J2iLi ®iiY I 
)e„ < iVmaxj(| dfjV and, as we have also assumed that V is a short range potential, 
the number of non-vanishing matrix elements d 2 jV is N(d + 1) where d is the number of 
neighbouring particles in the interaction range of the potential, thus (| d^d^V&V |) 
N(d+l)max i:j (\ frVdlVdW \) Sv . Finally, putting m = maxyfl d^dfjV^V |) s .„ 

55^ < maxid |) E „ | 2 _m(d+l) 



5j (.! 



< 



9f 



(9) 



which, in the limit iV — > oo, shows that the first derivative of the entropy is uniformly 
bounded by a finite constant. This first step proves that S^v) is continuous. 

The three further steps, concerning boundedness of the higher order derivatives, involve 
similar arguments to be applied to a number of terms which is rapidly increasing with the 
order of the derivative. But many of these terms can be grouped in the form of the variance 
or higher moments of certain quantities, thus allowing the use of a powerful technical trick 
to compute their iV-dependence. For example, using Eq.® in the expression for 8 2 Sn just 
below Eq.fJSJ), we get 



d 2 S 



N 



dv 



1 1 - 



< N 



(« 2 )s„- (a) 



+ N 



(10) 



where a = ||W|| A(1/||W||) and tp_ = V/||W||. Now, it is possible to think of the 
scalar function a as if it were a random variable, so that the first term in the r.h.s. of 
Eq. lJTUJ) would be its second moment. Such a possibility is related with the general validity 
of the Monte Carlo method to compute multiple integrals. In particular, since the H v are 
smooth, closed (V is non-singular), without critical points and re pres ent able as the union 
of suitable subsets of R^ -1 , the standard Monte Carlo method Jl|| is applicable to the 
computation of the averages (-)s„ which become sums of standard integrals in R^ -1 . This 
means that a random walk can be constructively defined on any £„, which conveniently 
samples the desired measure on the surface. Along such a random walk, usually called 
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Monte Carlo Markov Chain (MCMC), a and its powers behave as random variables whose 
"time" averages along the MCMC converge to the surface averages (-)sv Notice that the 
actual computation of these surface averages goes beyond our aim, in fact, we do not need the 
numerical values - but only the iV-dependences - of the upper bounds of the derivatives of 
the entropy. Therefore, all what we need is just knowing that in principle a suitable MCMC 
exists on each S v . Now, the function a is the integrand in square brackets in Eq.(J7J), where 
the second term vanishes at large N, as is clear from Eq.Q. Therefore, at increasingly 
large N, the approximate expression a = £) i=1 ^liY I II VV|| 2 tends to become exact, a is in 
the form of a sum function a = N^J^iLi^ °f terms a« = Nd^V/W VV^H 2 , of 0(1) in N, 
which, along a MCMC, behave as independent random variables with probability densities 
Ui(a,i) which we do not need to know explicitly. Then, after a classical ergodic theorem for 
sum functions, due to Khinchin ^(|, based on the Central Limit Theorem of probability 
theory, a is a gaussian-distributed random variable; as its variance decreases linearly with 
N, liniiv^oo N\(a 2 )z v - (a)|J = const < oo. 

Arguments similar to those above used for the first derivative of Sn lead to the result 
linijv^oo iV | V) ■ ip (a))sj = const < oo, which, together with what has been just found 
for the variance of a, proves the uniform boundedness also of the second derivative of Sn 
under the hypothesis of diffeomorphicity of the T, v . 

Similarly, but with an increasingly tedious work, we can treat the third and fourth deriva- 
tives of the entropy. In fact, despite the large number of terms contained in their expressions, 
they again belong only to two different categories: those terms which can be grouped in the 
form of higher moments of the function a, and whose iV-dependence is known after the 
above mentioned theorem due to Khinchin, and those terms whose iV-dependence can be 
found by means of the same kind of estimates given above for D^Sn- Eventually, after a 
lenghty but rather mechanical work, also the third and fourth derivatives of Sn are shown 
to be uniformly bounded as prescribed by Eq.pjl. Whence the proof of the Theorem. 

A few comments are in order. 

The converse of our Theorem is not true. There is not a one-to-one correspondence 
between phase transitions and topology changes, in fact, there are smooth, confining and 
finite-range potentials, like the one- dimensional XY model |7[ , with even a very large number 
of critical points, and thus many changes in the topology of the E„, but with no phase 
transition. Therefore, an open problem is that of sufficiency conditions, that is to determine 
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which kinds of topology changes can entail the appearance of a PT. Preliminary hints on 
this point are given by the analytic study of particular models [?], Il2| for which topology 
and thermodynamics are exactly computed. 

Finally, though at present our Theorem only applies to first and second order PTs and 
to those systems for which V is a good Morse function, it provides the grounding to an 
approach which can unify the mathematical description of very different kinds of PTs, like 
those "exotic" ones occurring in glasses or in the folding of polymers and proteins, for which 
the so-called energy landscape paradigm [3] is currently studied overlooking the link with 
Morse theory and topology. 

We thank H. van Beijeren, L. Casetti, E. Guadagnini, R. Lima, C. Liverani, A. Moro, P. 
Picco, M. Rasetti and G. Vezzosi for discussions and comments. 
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